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We investigate the dynamics of a two dimensional axial next nearest neighbour Ising 
(ANNNI) model following a quench to zero temperature. The Hamiltonian is given by _ff = 
— Jo ^.^j SijSi+ij — Ji ^^-^[SijSij+i — KSijSij+2]- For a < 1, the system does not reach the 
equilibrium ground state but slowly evolves to a metastable state. For k > 1, the system shows a 
behaviour similar to the two dimensional ferromagnetic Ising model in the sense that it freezes to a 
striped state with a finite probability. The persistence probability shows algebraic decay here with 
an exponent 6 = 0.235 ± 0.001 while the dynamical exponent of growth z — 2.08 ± 0.01. For k. = 1, 
the system belongs to a completely different dynamical class; it always evolves to the true ground 
state with the persistence and dynamical exponent having unique values. Much of the dynamical 
phenomena can be understood by studying the dynamics and distribution of the number of domains 
walls. We also compare the dynamical behaviour to that of a Ising model in which both the nearest 
and next nearest neighbour interactions are ferromagnetic. 

PACS numbers: 64.60.Ht, 75.60.Ch, 05.50.+q 



I. INTRODUCTION 



Dynamics of Ising models is a much studied phe- 
nomenon and has emerged as a rich field of present-day 
research. Models having identical static critical behav- 
ior may display different behavior when dynamic critical 
phenomena are considered [l|. An important dynamical 
feature commonly studied is the quenching phenomenon 
below the critical temperature. In a quenching process, 
the system has a disordered initial configuration corre- 
sponding to a high temperature and its temperature is 
suddenly dropped. This results in quite a few interest- 
ing phenomena like domain growth 0, persistence 
[i, 5, 6, 7„i etc. 

In one dimension, a zero temperature quench of the 
Ising model ultimately leads to the equilibrium configu- 
ration, i.e., all spins point up (or down). The average 
domain size D increases in time t as D{t) ^ t^/^, where 
z is the dynamical exponent associated with the growth. 
As the system coarsens, the magnetisation also grows in 
time as m{t) ~ t^/^^. In two or higher dimensions, how- 
ever, the system does not always reach equilibrium [H 
although these scaling relations still hold good. 

Apart from the domain growth phenomenon, another 
important dynamical behavior commonly studied is per- 
sistence. In Ising model, in a zero temperature quench, 
persistence is simply the probability that a spin has not 
flipped till time t and is given by P{t) ~ t"^ . 9 is called 
the persistence exponent and is unrelated to any other 
known static or dynamic exponents. 

Drastic changes in the dynamical behaviour of the Ising 
model in presence of a competing next nearest neighbor 
interaction have been observed earlier 0, [lo|, The 
one dimensional ANNNI (Axial next nearest neighbour 



Ising) model with L spins is described by the Hamiltonian 

L 

H = - J ^{S^Si+i - KSiSi+2). (1) 

i=l 

Here it was found that for k < 1, under a zero tempera- 
ture quench with single spin flip Glauber dynamics, the 
system does not reach its true ground state. (The ground 
state is ferromagnetic for k < 0.5, antiphase for k > 0.5, 
and highly degenerate at k = 0.5 [II])- On the contrary, 
after an initial short time, domain walls become fixed in 
number but remain mobile at all times thereby making 
the persistence probability go to zero in a stretched expo- 
nential manner. For k > 1 on the other hand, although 
the system reaches the ground state at long times, the dy- 
namical exponent and the persistence exponent are both 
different from those of the Ising model with only nearest 
neighbour interaction p^. 

The above observations and the additional fact that 
even in the two dimensional nearest neighbour Ising 
model, frozen-in striped states appear in a zero temper- 
ature quench 0], suggest that the two dimensional Ising 
model in presence of competing interactions could show 
novel dynamical behaviour. In the present work, we have 
introduced such an interaction (along one direction) in 
the two dimensional Ising model, thus making it equiv- 
alent to the ANNNI model in two dimensions precisely. 
The Hamiltonian for the two dimensional ANNNI model 
on a L X L lattice is given by 

L 

H = —Jo SijSi+ij~Ji [SijSijJi-i—KSijSij+2]- 

Henceforth, we will assume the competing interaction to 
be along the x (horizontal) direction, while in the y (ver- 
tical) direction, there is only ferromagnetic interaction. 
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Although the thermal phase diagram of the two dimen- 
sional ANNNI model is not known exactly, the ground 
state is known and simple. If one calculates the mag- 
netisation along the horizontal direction only, then for 
K < 0.5, there is ferromagnetic order and antiphase or- 
der for K > 0.5. Again, k = 0.5 is the fully frustrated 
point where the ground state is highly degenerate. On 
the other hand, there is always ferromagnetic order along 
the vertical direction. In Fig. [H we have shown the 
ground state spin configurations along the x direction 
for different values of k. 
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FIG. 2: The schematic pictures of configurations with flat 
interfaces separating domains of type I and II are shown: (a) 
when the interface lies parallel to y axis, we have nonzero fo^^ 
(= 2/L in this particular case) and (b) with interfaces parallel 
to the X axis we have nonzero foy (= 4/L here) 
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FIG. 1: The ground state (temperature T — 0) spin config- 
urations along the x direction are shown for different values 
of K. In the ferromagnetic phase, there is a two fold degen- 
eracy and in the antiphase the degeneracy is four fold. The 
ground state is infinitely degenerate at the fully frustrated 
point K — 0.5. 

In section II, we have given a list of the quantities cal- 
culated. In section III, we discuss the dynamic behaviour 
in detail. In order to compare the results with those 
of a model without competition, we have also studied 
the dynamical features of a two dimensional Ising model 
with ferromagnetic next nearest neighbour interaction, 
i.e., the model given by eq. ^ in which k < 0. These 
results are also presented in section III. Discussions and 
concluding statements are made in the last section. 



II. QUANTITIES CALCULATED 

We have estimated the following quantities in the 
present work: 

1. Persistence probability P{t): As already men- 
tioned, this is the probability that a spin does not 
flip till time t. 

In case the persistence probability shows a power 
law form, P{t) ~ , one can use the finite size 
scaling relation 



p{t,L)^r'f{Llt^/^) 



(3) 



For finite systems, the persistence probability satu- 
rates at a value at large times. Therefore, for 
X << 1 , /(x) ^ x~°' with a = z9. For large x, 
f{x) is a constant. 

It has been shown that the exponent a is related to 
the fractal dimension of the fractal formed by the 
persistent spins [l3j . Here we obtain an estimate of 
a using the above analysis. 



Number of domain walls : Taking a single strip 
of L spins at a time, one can calculate the number 
of domain walls for each strip and determine the 
average. In the L x L lattice, we consider the frac- 
tion fo — No/L and study the behaviour of fu as 
a function of time. One can take strips along both 
the X and y directions (see Fig. [5] where the calcu- 
lation of fu in simple cases has been illustrated). 
As the system is anisotropic, it is expected that the 
two measures, fo^ along the x direction and foy 
along the y direction, will show different dynamical 
behaviour in general. The domain size D increases 
as t^^^ as already mentioned and it has been ob- 
served earlier that the dynamic exponent occurring 
in coarsening dynamics is the same as that occur- 
ring in the finite size scaling of P{t) (eq. ([3])) p^ . 
Although we do not calculate the domain sizes, the 
average number of domain walls per strip is shown 
to follow a dynamics given by the same exponent 
z, at least for k > 1. 



3. Distribution Pifo) (or P{Nd)) of the fraction (or 
number) of domain walls at steady state: this is 
also done for both x and y directions. 



4. Distribution P{m) of the total magnetisation at 
steady state for k < only. 

We have taken lattices of size L x L with L — 
40, 100, 200 and 300 to study the persistence be- 
haviour and dynamics of the domain walls of the 
system and averaging over at least 50 configura- 
tions for each size have been made. For estimating 
the distribution Nu we have averaged over much 
larger number of configurations (typically 4000) 
and restricted to system sizes 40 x 40, 60 x 60, 
80 X 80 and 100 x 100. Periodic boundary con- 
dition has been used in both x and y directions. 
Jo = Ji = I has been used in the numerical simu- 
lations. 
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III. DETAILED DYNAMICAL BEHAVIOUR 

Before going in to the details of the dynamical 
behaviour let us discuss the stability of simple config- 
urations or structures of spins which will help us in 
appreciating the fact that the dynamical behaviour is 
strongly dependent on k. 



A. Stability of simple structures 

An important question that arises in dynamics is the 
stability of spin configurations - it may happen that con- 
figurations which do not correspond to global minimum 
of energy still remain stable dynamically. This has been 
termed "dynamic frustration" earlier. A known ex- 
ample is of course a striped state occurring in the two or 
higher dimensional Ising models which is stable but not 
a configuration which has minimum energy. 

In ANNNI model, the stability of the configurations 
depend very much on the value of k. It has been previ- 
ously analysed for the one dimensional ANNNI model 
that K = 1 is a special point above and below which the 
dynamical behaviour changes completely because of the 
stability of certain structures in the system. 

Let us consider the simple configuration of a single up 
spin in a sea of down spins. Obviously, it will be unstable 
as long as K < 2. For k > 2, although this spin is stable, 
all the neighbouring spins are unstable. However, for 
K < 2, only the up spin is unstable and the dynamics 
will stop once it flips. When k = 2 the spin may or may 
not flip, i.e., the dynamics is stochastic. 

Next we consider a domain of two up spins in a sea of 
down spin. These two may be oriented either along hori- 
zontal or vertical direction. These spins will be stable for 
K > 1 only while all the neighouring spins are unstable. 
For K < 1, all spins except the up spins are stable. When 
K = 1, the dynamics is again stochastic. 

A two by two structure of up spins in a sea of down 
spins on the other hand will be stable for any value of 
K > 0. But the neighbouring spins along the vertical 
direction will be unstable for k > 1. This shows that 
for K < 1, one can expect that the dynamics will affect 
the minimum number of spin and therefore the dynamics 
will be slowest here. A picture of the structures described 
above are shown in Fig[3l 

One can take more complicated structures but the 
analysis of these simple ones is sufficient to expect that 
there will be different dynamical behaviour in the regions 
K < 1,K = 1,K > 1,K = 2 and k > 2. However, we 
find that as far as persistence behaviour is concerned, 
there are only three regions with different behaviour: 
K < 1, K — 1 and K > 1. On the other hand, when the 
distribution of the number of domain walls in the steady 
state is considered, the three regions 1 < k < 2, k — 2 
and K > 2 have clearly distinct behaviour. 
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FIG. 3: Analysis of stability of simple structures: (a) single 
up spin in sea of down spins; here for k < 2 all the spins 
except the up spin is stable (b) two up spins in a sea of down 
spins, all spins except the two up spins are stable for k < 1 
(c) a two by two structure of up spins - here all the spins are 
stable for k < 1 while neighbouring spins are not (see text for 
details) . 



B. < K < 1 



We find that as in [l^, in the region < k < 1, the 
system has identical dynamical behaviour for all k. Also, 
like the one dimensional case, here the system does not go 
to its equilibrium ground state. However, the dynamics 
continues for a long time, albeit very slowly for reasons 
mentioned above. In Figs. |4]-[7l we show the snapshots of 
the system at different times for a typical quench to zero 
temperature. As already mentioned, here domains of size 
one and two will vanish very fast and certain structures, 
the smallest of which is a two by two domain of up/down 
spins in a sea of oppositely oriented spins can survive till 
very long times. These structures we call quasi-frozen 
as the spins inside these structures (together with the 
neighbourhood spins) are locally stable; they can be dis- 
turbed only when the effect of a spin ffip occurring at a 
distance propagates to its vicinity which usually takes a 
long time. 

The pictures at the later stages also show that the sys- 
tem tends to attain a configuration in which the domains 
have straight vertical edges, it can be easily checked that 
structures with kinks are not stable. We find a tendency 
to form strips of width two ( "ladders" ) along the vertical 
direction - this is due to the second neighbour interac- 
tion - however, these strips do not span the entire lattice 
in general. The domain structure is obviously not sym- 
metric, e.g., ladders along the horizontal direction will 
not form stable structures. The dynamics stops once the 
entire lattice is spanned by only ladders of height M < L. 



FIG. 4: Snap shot of a 40 x 40 system at time t — 10 for 
K < 1 A few simplest quasi frozen structures are highlighted. 



FIG. 6; Same as Fig. g] with t = 500. One of the two by 
two structures has melted while another one has formed. The 
ladder like structures which have formed are perfectly stable. 
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FIG. 5: Same as Fig. Hwith t = 100. 



The persistence probability for k < 1 shows a very 
slow decay with time which can be approximated by -j^^^ 
for an appreciable range of time. At later times, it ap- 
proaches a saturation value in an even slower manner. 
The slow dynamics of the system accounts for this slow 
decay. 

The fraction of domain walls fjj^ and Jd^ along the 
X direction and y directions show remarkable difference 
as functions of time. While that in the x direction satu- 
rates quite fast, in the y direction, it shows a gradual de- 
cay till very long times (see Fig. [5]). This indicates that 
the dynamics essentially keeps the number of domains 
unchanged along x direction while that in the other di- 
rection changes slowly in time. The behaviour of fo^ is 
similar to what happens in one dimension. In fact, the 
average number of domain walls at large times is 

also very close to that obtained for the ANNNI chain, it 
is about 0.27L. However, in contrast to the one dimen- 
sional case where the domain walls remain mobile, here 
the mobility of the domain walls are impeded by the pres- 
ence of the ferromagnetic interaction along the vertical 
direction causing a kind of pinning of the domain walls. 

The distribution of the fraction of domain walls in the 
steady state shown in Fig. [9] also reveals some important 
features. The distribution for fo^ and fu are both quite 
narrow with the most probable values being fu^ ~ 0.27 
and foy — 0.04 (these values are very close to the av- 
erage values). With the increase in system size, the dis- 
tributions tend to become narrower, indicating that they 
approach a delta function like behaviour in the thermo- 
dynamic limit. 




FIG. 7: Same as Fig. g] with t = 75000. This snapshot 
is taken after a very long time to show that the system has 
undergone nominal changes compared to the length of the 
time interval. The whole configuration now consists of ladders 
and the dynamics stops once the system reaches such a state. 

C. K> 1 

It was already observed that k = 1 is the value at which 
the dynamical behaviour of the ANNNI model changes 
drastically in one dimension. In two dimensions, this is 
also true, however, wc find that the additional ferromag- 
netic interaction along the vertical direction is able to 
affect the dynamics to a large extent. Again, similar to 
the one dimensional case, we have different dynamical 
behaviour for k = 1 and k > 1. In this subsection we 
discuss the behaviour for k > 1 while the k — 1 case is 
discussed in the next subsection. 

The persistence probability follows a power law decay 
with 6 = 0.235 ± 0.001 for all k > 1, while the finite 
size scaling analysis made according to ([3|) suggests a z 
value 2.08 ± 0.01. This is checked for different values of 
K {k = 1.3,1.5,2.0,20,100) and the values of 9 and z 
have negligible variations with n which do not show any 
systematics. Hence we conclude that the exponents are 
independent of k for k > 1. A typical behaviour of the 
raw data as well as the data collapse is shown in Fig. [TOl 

The dynamics of the average fraction of domain walls 
along the horizontal direction, fo^ again shows a fast 
saturation while that in the y direction has a power law 
decay with an exponent ~ 0.48 (Fig. [TT|) . This exponent 
is also independent of k. As mentioned in section II, 
we find that there is a good agreement of the value of 
this exponent with that of 1/z obtained from the finite 
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FIG. 8: Persistence P{t) and average number of domain walls 
per site, fo are shown for k < 1. 
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FIG. 9: Steady state distributions of fraction of domain walls 
at K < 1 for different system sizes. The distributions become 
narrower as the system size is increased. 



Interfaces which occur parallel to the y axis, separating 
two regions of antiphase and keeping the ferromagnetic 
ordering along the vertical direction intact, are extremely 
rare, the probability vanishing for larger sizes. Quan- 
titatively this means we should get Ju^ = 0.5 at long 
times which is confirmed by the data (Fig. Ilip . Hence in 
the following our discussions on striped state will always 
imply flat horizontal interfaces, i.e., antiphase ordering 
along each horizontal row but the ordering can be of dif- 
ferent types (e.g., a-l-H V ^ •■• type and a 

V ^ !- + ••• type, which one can call a 'shifted' 

antiphase ordering with respect to the first type) . 

It is of interest to investigate whether these striped 
states survive in the infinite systems. To study this, we 
consider the distribution of the number of domain walls 
rather than the fraction for different system sizes. The 
probability that there are no domain walls, or a perfect 
ferromagnetic phase along the vertical direction, turns 
out to be weakly dependent on the system sizes but hav- 
ing different values for different ranges of values of k. For 
1 < K < 2, it is ~ 0.632, for k = 2.0, it is ~ 0.544 while 
for any higher value of k, this probability is about 0.445. 
Thus it increases for k although not in a continuous man- 
ner and like the two dimensional case, we find that there 
is indeed a finite probability to get a striped state. 

While we look at the full distribution of the number of 
domain walls at steady state (Fig. [T5)l . we find that there 
are dominant peaks at No = (corresponding to the 
unstriped state) and at = 2 (which means there are 
two interfaces). However, we find that the distribution 
shows that there could be odd values of as well. 
This is because the antiphase has a four fold degeneracy 
and the and a 'shifted' ordering can occur in several ways 
such that odd values of are possible. In any case, 
the number of interfaces never exceeds iVn = 6 for the 
system sizes considered. 



size scaling behaviour of P{i) implying that the average 
domain size D is inversely proportional ■ This is quite 
remarkable, as the fraction of domain walls calculated in 
this manner is not exactly equivalent to the inverse of 
domain sizes in a two dimensional lattice; the fact that 
/o^ remains constant may be the reason behind the good 
agreement (essentially the two dimensional behaviour is 
getting captured along the dimension where the number 
of domain walls show significant change in time) . 

Although the persistence and dynamic exponents are n 
independent, we find that the distribution of the number 
of domain walls has some nontrivial n, dependence. 

Though the system, for all k > 1, evolves to a state 
with antiphase order along the horizontal direction, the 
ferromagnetic order along vertical chains is in some cases 
separated by one or more domain walls. A typical snap- 
shot is shown in Fig. [T^] displaying that one essentially 
gets a striped state here like in the two dimensional Ising 
model. 
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FIG. 10: The collapse of scaled persistence data versus scaled 
time using Q — 0.235 and z — 2.08 is shown for different 
system sizes for k > 1. Inset shows the unsealed data. 
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FIG. 11: Decay of the fraction of domain walls with time at 
K > 1 are shown along horizontal and vertical directions. The 
dashed line has slope equal to 0.48. 
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FIG. 13: Normalised steady state distributions of number 
of domain walls for different k > 1 show that striped states 
occur with higher probability as k increases. The lines are 
guides to the eye. 
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FIG. 12: A typical snapshot of a steady state configuration 
for K > 1 with fiat horizontal interfaces separating two regions 
of antiphase ordering (see text). 




FIG. 14: The collapse of scaled persistence data versus scaled 
time using 6 — 0.263 and z = 1.84 is shown for different 
system sizes at k = 1. Inset shows the unsealed data. 



D. K = l 

Here we find that the persistence probability follows a 
power law decay with 6 = 0.263 ± 0.001. The finite size 
scaling analysis suggests a z value 1.84 ± 0.01 (Fig. [T4|) . 

We have again studied the dynamics of fjj^ and fjj^ ; 
the former shows a fast saturation at 0.5 while the latter 
shows a rapid decay to zero after an initial power law 
behaviour with an exponent « 0.515 (Fig. [T5|) . This 
value, unlike in the case k > 1, does not show very good 
agreement with 1/z obtained from the finite size scal- 
ing analysis. We will get back to this point in the next 
section. 

The results for fu^ and foy imply that the system 
reaches a perfect antiphase configuration as there are no 
interfaces left in the system with fu^ — 0.5 and fjj = 
at later times. 



E. K< 0.0 

In order to make a comparison with the purely ferro- 
magnetic case, we have also studied the Hamiltonian ([2|) 
with negative values of k which essentially corresponds 
to the two dimensional Ising model with anisotropic next 
nearest neighbour ferromagnetic interaction. 

K — corresponds to the pure two dimensional Ising 
model for which the numerically calculated value of ~ 
0.22 is verified. We find a new result when k is allowed to 
assume negative values, the persistence exponent 9 has a 
value ~ 0.20 for |k| > 1 while for < |k| < 1, the value of 
6 has an apparent dependence on k, varying between 0.22 
to 0.20. However, it is difficult to numerically confirm 
the nature of the dependence in such a range and we 
have refrained from doing it. At least for |k| >> 1, the 
persistence exponent is definitely different from that of 
at K = 0. The growth exponent z however, appears to 
be constant and ~ 2.0 for all values of k < 0. A data 
collapse for large negative k is shown in Fig. [12] using 
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FIG. 15: Decay of the fraction of domain walls with time at 
n — 1 are shown along horizontal and vertical directions. The 
dashed line has slope equal to 0.515. 
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FIG. 17: Decay of the fraction of domain walls with time 
at «: < are shown along horizontal (fo^), shown by dotted 
lines) and vertical (/d„), shown by solid lines) directions. 
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FIG. 16: The collapse of scaled persistence data versus scaled 
time using 9 — 0.20 and z — 2.0 is shown for different system 
sizes for «: < — 1. Inset shows the unsealed data. 

The effect of the anisotropy shows up clearly in the 
behaviour of Jd^ and fn as functions of time (Fig. \T7\ . 
For K = 0, they have identical behaviour, both reaching 
a finite saturation value showing that there may be inter- 
faces generated in either of the directions (corresponding 
to the striped states which are known to occur here) . As 
the absolute value of k is increased, fo^. shows a fast de- 
cay to zero while foy attains a constant value. The sat- 
uration value attained by foy increases markedly with 
|k| while for fu^ the decay to zero becomes faster. One 
can conduct a stability analysis for striped states to show 
that such states become unstable when the interfaces are 
vertical and k increases beyond 1, leading to the result 
Id^ ^ 0. 

Extracting the z value from the variations of fjj^ or 
fo is not very simple here as the quantities do not show 
smooth power law behaviour over a sufficient interval of 
time. 

The fact that and/or fu^ reach a finite saturation 



value indicates that striped states occur here as well. The 
behaviour of fo^ and fu suggests that in contrast to the 
isotropic case where interfaces can appear either horizon- 
tally or vertically, here the interfaces appear dominantly 
along the x direction as k is increased. Thus the nor- 
malised distribution of the number of domain walls along 
y is shown in Fig. [THl We find that as k is increased in 
magnitude, more and more interfaces appear. However, 
the number of interfaces is always even consistent with 
the fact that interfaces occur between ferromagnetic do- 
mains of all up and all down spins. 
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FIG. 18: Normalised steady state distributions of number 
of domain walls for different k < show that striped states 
occur with higher probability as |k| increases. The lines are 
guides to the eye. 

Lastly in this section, we discuss the behaviour of the 
magnetisation which is the order parameter in a fer- 
romagnetic system. As striped states are formed, the 
magnetisation will assume values less than unity. The 
probability of configurations with magnetisation equal to 
unity shows a stepped behaviour, with values changing at 
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= 1 and 2 and assuming constant values at 1 < |k| < 2 
and above |k| = 2 (Fig. [H]). 
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FIG. 19: Probability that the magnetisation takes a steady 
state value equal to unity is shown against k when ac < 0. 



IV. DISCUSSIONS AND CONCLUSIONS 

We have investigated some dynamical features of the 
ANNNI model in two dimensions following a quench to 
zero temperature. We have obtained the results that the 
dynamics is very much dependent on the value of k, the 
ratio of the antiferromagnetic interaction to the ferro- 
magnetic interaction along one direction. This is similar 
to the dynamics of the one dimensional model studied 
earlier, but here we have more intricate features, e.g., 
that of the occurrence of quasi frozen-in structures for 
K < 1 where the persistence probability shows a very slow 
decay with time. Persistence probability is algebraic for 
K > 1, but exactly at k = 1, the exponents 9 and z are 
different from those at k > 1. The exponents for k > 1 
are in fact very close to those of the two dimensional Ising 
model with nearest neighbour ferromagnetic interaction. 
(This was not at all true for the one dimensional ANNNI 
chain, where the persistence exponent at k > 1 was found 
to be appreciably different from that of the one dimen- 
sional Ising chain with nearest neighbour ferromagnetic 
interaction.) This shows that the ferromagnetic inter- 
action along the vertical direction is able to negate the 
effect of the antiferromagnetic interaction to a great ex- 
tent. This is apparently a counter intuitive phenomenon, 
K — and K > 1 having very similar dynamic behaviour 
while in the intermediate values, the dynamics is qualita- 
tively and quantitatively different. At far as dynamics is 
concerned, the ANNNI model in two dimensions cannot 
be therefore treated perturbatively. 

Although the values of 9 and z are individually quite 
close for k — and k > 1, the product z9 — a are quite 
different. For k = 0, a ~ 0.44 while for k > 1, it is 
0.486 ± 0.002. This shows that the spatial correlations 
of the persistent spins are quite different for the two and 



one can safely say that the dynamical class for k = 
and K > 1 are not the same, k = 1 is the special point 
where the dynamic behaviour changes radically. Here 
there appears to be some ambiguity regarding the value 
of z; estimating a from the finite size scaling analysis 
gives a « 0.484 ± 0.005 while using the z value from the 
domain dynamics, the estimate is approximately equal 
to 0.51. However, the dynamics of the domain sizes may 
not be very accurately reflected by the dynamics of foy 
in which case a w 0.48 is a more reliable result. Thus 
we find that although the values of 9 and z are quite 
different for k = 1 and k > 1, the a values are close. 

We would like to add here that when there is a power 
law decay of a quantity related to the domain dynamics, 
it is highly unlikely that it will be accompanied by an 
exponent which is different from the growth exponent. 
Thus, even though we get slightly different values of z 
for K = I from the two analyses, it is more likely that 
this is an artifact of the numerical simulations. 

Another feature present in the two dimensional Ising 
model is the finite probability with which it ends up in a 
striped state. The same happens for k > 1, but here the 
probabilities are quite different and also dependent on k. 
We find that there is a significant role of the point k — 2 
here as this probability has different values at k — 2, 
K > 2 and k < 2. 

Comparison of the ANNNI dynamics with that of the 
ferromagnetic anisotropic Ising model shows some inter- 
esting features. In the latter, one gets a new value of 
persistence exponent for k < — 1 while in the former a 
new value is obtained for k > l.The new values (except 
for K = 1) are in fact very close to that of the two di- 
mensional Ising model, but simulations done for identi- 
cal system sizes averaged over the same number of initial 
configurations are able to confirm the difference. The 
qualitative behaviour of the domain dynamics is again 
strongly k dependent when n is negative. Another point 
to note is that the probability that the system evolves to a 
pure state is k dependent in both the ANNNI model and 
the Ising model. In both cases in fact, this probability 
decreases in a step like manner with increasing magni- 
tude of K. We also find the interesting result that while 
the distribution of the number of domain walls can have 
non-zero values at odd values of No in the ANNNI model 
because of the four fold degeneracy of the antiphase, for 
the Ising model, odd values of Njj are not permissible as 
the ferromagnetic phase is two fold degenerate. 

Finally we comment on the fact that although the dy- 
namical behaviour, as far as domains are concerned, re- 
flects the inherent anisotropy of the system (in both the 
ferromagnetic and antiferromagnetic models), the persis- 
tence probability is unaffected by it. In order to verify 
this, we estimated P{t) along an isolated chain of spins 
along X and y directions separately and found that the 
two estimates gave identical results for all values of k. 

In conclusion, it is found that except for the region 
< |k| < 1, the dynamical behaviour of the Hamiltonian 
^ is remarkably similar for negative and positive k; the 
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persistence and growth exponents get only marginally 
affected compared to the values of the two dimensional 
Ising case (k = 0) and the domain distributions have 
similar nature. However, the region < k < 1 is extraor- 
dinary, where algebraic decay of persistence is absent. 
There is dynamic frustration as the system gets locked 
in a metastable state consisting of ladder-like domains 
and the dynamics is very slow because of the presence 
of quasi-frozen structures. There is in fact dynamic frus- 
tration at other k values also in the sense that except 
for K = 1, the system has a tendency to get locked in 
a "striped state". However, even in that case, the al- 
gebraic decay of the persistence probability is observed. 
Thus algebraic decay of persistence probability seems to 
be valid only when the metastable state is a striped state. 
Although there is no dynamic frustration at k = 1 in the 
sense that it always evolves to a state with perfect an- 
tiphase structure, it happens to be a very special point 
where the persistence exponent and growth exponents 
are unique and appreciably different from those of the 
K = case. 

In this paper, the behaviour of the two dimensional 
ANNNI model under a zero temperature has been dis- 



cussed; the dynamics at finite temperature can be in fact 
quite different. At finite temperatures, the spin flipping 
probabilities are stochastic, and the dynamical frustra- 
tion may be overcome by the thermal fluctuautions. It 
has been observed earlier [ij] that in a thermal anneal- 
ing scheme of the one dimensional ANNNI model, the 
K = 0.5 point becomes significant. A similar effect can 
occur for the two dimensional case as well. The defini- 
tion of persistence being quite different at finite temper- 
atures [1^ , it is also not easy to guess its behaviour (for 
either the one or two dimensional model) simply from 
the results of the zero temperature quench. Indeed, the 
ANNNI model under a finite temperature quench is an 
open problem which could be addressed in the future. 
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